Gravitational Lens Time Delays in CDM 
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ABSTRACT 

In standard CDM halo models, the time delay of a gravitational lens is determined 
by the cold baryon mass fraction, /& = £lb,coid/^o, of the visible galaxy relative 
to the overall halo. The observed time delays in PG1115+080, SBS1520+530, 
B1600+434 and HE2149-2745 give Hubble constants consistent with the HST Key 
Project value of Hq = 72 ± 8 km/s Mpc only if /& <; 0.2 (1-sided 68% confidence), 
which is larger than the upper bound of fb,max = ^fe/^o = 0.15 ± 0.05 estimated 
from the CMB. If all available baryons cool and /& = fb, ma x then the time delays 
imply Hq = 65 ± 6 km/s Mpc (95% confidence). If local inventories of cold baryons, 
fb — 0.013//&70, are correct, then Hq = 52 ± 6 km/s Mpc and the halo parameters 
closely match isothermal mass models. Isothermal models are also consistent with 
strong and weak lens studies, stellar dynamics and X-ray observations on these scales, 
while significantly more centrally concentrated models are not. There is a a conflict 
between gravitational lens time delays, the local distance scale and standard CDM 
halo models. 



Subject headings: cosmology: gravitational lensing; cosmology: Hubble constant; dark 
matter 



Introduction 



Kochanek (2002a) found that it was difficult to reconcile the time delays measured for 5 
simple, well-observed gravitational lenses with the local distance scale given our expectation that 
galaxies have massive, extended dark matter halos. If the lens galaxies had constant mass-to-light 
(M/L) ratios we found Hq = 71 =L 6 km/s Mpc, which is consistent with the local estimate 
of Hq = 72 ± 8 km/s Mpc by the HST Key Project (Freedman et al. 2001 ). However, if the 



lenses had isothermal mass distributions (flat rotation curves), we found Hq = 48_ 4 km/s Mpc, 
which is grossly inconsistent with the HST Key Project. While the time delay lenses cannot 
distinguish between these two limiting mass distributions, models of other lenses (e.g. Munoz, 
Kochanek & Keeton 2001 ), stellar dynamical measurements (e.g. Rix et al. 1997| , Romanowsky 
& Kochanek 199£, Gerhard et al. [2001 , Treu &; Koopmans |2002| ), weak lensing (e.g. Guzik Sz 
Seljak 2002) and X-ray (e.g. Fabbiano 198S, Lowenstein & White |1999 ) measurements all suggest 
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that the isothermal mass distributions are correct. In this study we will show that standard cold 
dark matter (CDM) halo models closely resemble the isothermal models on these scales, which 
implies there is a conflict between the local distance scale, gravitational lens time delays and CDM 
halo models. 



While the Kochanek (2002a) results provided evidence for a real conflict given the considerable 
observational evidence that lens galaxies must have extended, massive dark matter halos, the link 
to a problem with CDM halo models was qualitative because the study lacked a quantitative, 
theoretical prediction for the time delays expected from CDM halos. One barrier to making such 
predictions was that we lacked a clear understanding of which features of lens mass distributions 
control time delays. While global degeneracies due to the addition of constant mass density sheets 
(e.g. Falco, Gorenstein & Shapiro 1985 , Gorenstein, Falco & Shapiro 1988] , Saha 2000 ) and a 
correlation between more compact mass distributions and longer time delays (e.g. Schechter 2000| . 



Witt, Mao & Keeton 2000) were well known, it was unclear which properties of a halo model 
had to be accurately computed in order to make robust predictions. Kochanek ( |2002b ) combined 
analytic results with comparisons to the numerical models by Kochanek ( |2002a ) to show that 
the surface density in the annulus between the images used to measure the delay was the most 
important physical property of the lens galaxy for determining the time delay. The interior mass 
is implicit in the astrometry of the images and the lens galaxy, and the angular structure is either 
unimportant or strongly constrained by the astrometry. As a result, the Hubble constant expected 
for a simple lens is related to the surface density by Hq = A(l — (k)) + B(r] — 1)(k), where (k) is the 
average surface density in the annulus between the images (in units of the critical density), with 
a modest correction \B\ <, A/ 10 due to the logarithmic slope rj of the surface density distribution 
within the annulus (k oc R}~~' n ). The coefficients A and B depend only on the image positions and 
the measured time delay. These simple semi-analytic scaling laws reproduce full numerical models 
to accuracies of better than 5%. 

We can now calculate the expected properties of gravitational lens time delays for CDM halo 
models. In §2 we outline our model for the halos, which are based on the CDM lens models from 
Keeton ( 2001 ). The models consist of a Hernquist ( fl990| ) model for the luminous early- type lens 
galaxy embedded in an NFW (Navarro, Frenk & White |1996 ) halo normalized using the parameter 
estimates of Bullock et al. (|2001[ ). We considered both unmodified NFW halos and adiabatically 
compressed (Blumenthal et al. 1986| ) halos. We summarize the mathematical details of the 
model in §2. In §3 we apply it to the four simple time delay lenses PG1115+080, SBS1520+530, 
B1600+434 and HE2149-2745, to show that the values of («) and rj that determine the Hubble 
constant given the measured time delays are in turn determined by a single parameter, the cold 
baryonic mass fraction, = £lb,coid/^o, °f the luminous galaxy compared to the halo. Since the 
baryon fraction is bounded by local estimates from observed baryonic populations and the global 
baryon fraction estimated either in clusters or from the CMB, we can set firm bounds for the 
range of Hq consistent with CDM halo models. As we discuss in §4, this leads to a new element 



of the so-called "dark matter crisis" (e.g. Moore |200l[) , because the CDM halo models combined 
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with the measured time delays require lower Hubble constants than are consistent with the Key 
Project estimates based on the local distance scale. An Appendix briefly discusses the effects of 
tidal truncation on lens galaxy halos. 



2. The CDM Lensing Model 

Our model for lensing by galaxies with standard CDM halos follows that of Keeton ( |2001| ). 
These models are very similar to the lens models used by Kochanek & White (2001) or those 
used in semi-analytic models for the structure of galaxies (e.g. Mo, Mao & White 1998, Cole et 
al. pOOOj , Gonzalez et al. 2000 and references therein), but focus on early- type rather than late- type 
galaxies. We model the visible lens galaxy with a Hernquist ( 199C| ) profile, 



Ph{t) 



r H 



2tt r{r + r#) 3 ' 



(1) 



with the scale length of rn = 0.551i? e that matches the Hernquist profile to a de Vaucouleurs 
profile with effective radius R e . The enclosed mass of the Hernquist model is 



M H (< r) = M H - 



(r + r H ) 2 ' 

We model the initial dark matter halo using an NFW (Navarro, Frenk & White |l99"6| ) profile, 



Pn(t) 



M„ 



47r/(c) r(r + r s ) s 



(2) 



(3) 



normalized by the mass M v i r inside the radius r v i r with concentration c = r v i r /r s and 
/(c) = ln(l + c) — c/(l + c). The enclosed mass is 

M N (<r)=M vir f(r/r s )/f(c). 



(4) 



:i~ — 1| 1 / 2 tann 1 \x 2 — l] 1 / 2 , where the tann 1 — > tan 1 (tanh 1 ) for x > 1 {x < 1), then 



We need the properties of these potentials in projection for the lensing calculations. If we define 

F(x) = \x 2 

the projected surface density of the Hernquist profile is 

M H (2 + x 2 )T{x) - 3 



T, H (x = R/r H ) 
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and the mass inside cylindrical radius R is 

M(< x = R/rjj) = M H 



{X 2 - l) 2 
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(Hernquist 199C, Keeton 2001). Similarly, the projected properties of the NFW model are 

M vir 1-T{x) M vir , 



T,n(x = R/r s 



27rr 2 /(c) x 2 - 1 



-Etv(x) 



(5) 
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and 



M N (< x = R/r s 



M„ 



ln| + JT( x ) 



M vir M N (< x) 



(8) 



for the projected surface density and mass respectively (Bartelmann |1996 ), The lensing properties 
of the halo depend on the surface density measured in units of the critical surface density, 
K = S/S C5 where the critical surface density, S c = c 2 Dos I '^kGDolDls > is a simple function of 
the angular diameter distances between the observer, the lens and the source (e.g. Schneider, 



Ehlers & Falco |1992|) 



We need only two observational parameters to normalize the models. The first is the 
(intermediate axis) effective radius R e of the lens galaxy, which determines the density distribution 
of the visible baryons (stars). The second is the critical radius R c of the lens, determined either 
by detailed models or simply by the average distance of the images from the lens center, which 
determines the mass scale. The halo properties depend on the halo concentration, c, virial mass 
and radius, M v i T and r„j r , and the mass fraction, /j, represented by the luminous galaxy. We 
adopt the halo parameter estimates from Bullock et al. ( |2001| ) and an Qq = 0.3 flat cosmological 
model. The virial mass and radius are related by a definition, 



4tt 
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where p u (z) = 3HqQq(1 + z) 3 /8ttG is the mean matter density at the lens redshift and the virial 
overdensity is A v i r ~ (18-7T 2 + 82x — 39x 2 )/£l(z) where x = ft(z) — 1. Simulations find that the 
break radius and the virial radius are statistically correlated through the concentration parameter 
c = r v i r /r s . The average concentration depends on the halo mass M v ^ r and the redshift 
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-0.14 



(10) 



but individual halos have a log-normal dispersion about the average concentration of approximately 
cr c = 0.18 (base 10). Finally, a fraction of the mass cools to form the visible galaxy, so the mass 
of the Hernquist profile is f\,M v i r and the remaining mass of the halo is (1 — fb)M v i r . 

We use the geometry of the lens to determine the overall mass of the halo. The total projected 

= 7T.R 2 E C (e.g. Schneider et al. |1992|) , 



mass inside the critical radius of a spherical lens is M(< R, 
or 



vri? 2 i; r = M v 



f b M H (Rc/r H ) + (1 - f b )M N {Rc/r s 



(11) 



for a model without adiabatic compression. We combine this equation with the relation between 
the virial mass and radius (Eqn. ||) to determine M v i r and r v i r given the cold baryon fraction /{, 
and the concentration c. The surface density, in units of the critical density, is then 



k(R) 



h-£ H (R/r H ) + - 



fb 



X N (R/r s 



(12) 
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from which we can estimate both the average surface density (re) and the logarithmic slope of the 
surface density rj (assuming k oc i? 1_r? locally) in the annulus between the lensed images. 

These models probably underestimate the central density of the dark matter because 
they neglect the compression of the dark matter by the cooling of the baryons. We use the 



adiabatic compression model (Blumenthal et al. 1986) to compute the changes in the dark matter 



distribution. For a spherical system with particles on circular orbits, slow changes in the radial 
mass distribution preserve the angular momentum of the particles. For initial and final radii of fj 
and rj and mass distributions Mj(< r) and Mf(< r), the angular momentum of a circular orbit is 
conserved if rjMj(< rj) = rfMf(< r/). The initial mass distribution is simply that of the NFW 
halo, Mj(< r) = Mn(< r). The final mass distribution is 

M f(< r f) = fbM H (< r f ) + (1 - f b )M A (< Tf ) = f b M H (< r f ) + (1 - f h )M N {< n) (13) 

where the final halo mass distribution M A {< rf) = M v i r M A (< rf is the same as that of the 
original halo at the initial radius Mjy(< rj). For the compression of a NFW halo by Hernquist 
galaxy, the solution for rj(rj) is analytic (a cubic equation, see Keeton 2001 ). Jesseit et al. ( [2002; ) 



show that simple Blumenthal et al. (1986) model accurately estimates the changes in the mass 
distribution. The compressed halo density profile must then be numerically projected to determine 
the projected mass 

M A {< R) = M vir M A (< R) = M vir / d(j)COS 4>M A {< r = Rsec 4>) (14) 

Jo 

and surface density 



E A {R) = — —Z A (R) = —— d(p——{<r = Rsec(f>), (15) 
ri ZttK Jn ar 



and then we must solve Eqns. (|ll]) and fli~2| ) for M v i r and r v i r with the NFW profile replaced by 
the adiabatically compressed profile. For one lens, B1600+434, we include an exponential disk 
as well as a bulge in the baryonic model, fixing the mass ratio to that inferred from the constant 



M/L mass models of Kochanek ( 2002a ). Once we have determined the surface density, we can 
easily calculate the mean surface density (k) and the logarithmic slope r] of the surface density in 
the annulus between the images for which time delays have been measured. 

In addition to estimates of the Hubble constant, we can also calculate the weak lensing 
aperture mass, AE(i?) = (£(< R)) — S(i?), from the surface density profile. Because weak lensing 
measurements are made on much larger scales (R J> 50h^ 1 kpc) than the critical radius of the 
lens (R ~ 5/1" 1 kpc), they provide a nearly ideal additional constraint on the extent and mass 
of the dark matter halo. Compact, low-mass, constant M/L models will produce smaller weak 
lensing signals than extended, high-mass, dark matter dominated models. While there are no 
weak lensing measurements for time delay lenses, we can estimate the expected signal from the 



measurements by the Sloan Digital Sky Survey (SDSS, McKay et al. |2001| , Guzik & Seljak |2002|) . 
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Fig. 1. — The average surface density (k) for adiabatically compressed models of HE2149-2745 as 
a function of the halo concentration c = r v i r /r s and the mass fraction in the cold baryons of the 
visible galaxy = Mn/M v i r . The solid lines are contours of (k) spaced at intervals of A(k) = 0.1. 
An isothermal model has (k) = 0.5 (the heavy solid contour). Constant M/L models lie on the top 
edge where /& = 1. In the Bullock et al. ( 2001D simulations of halo formation, halos have a limited 
range of concentrations for a given virial mass. The heavy dashed contour shows the most likely 
halo concentration, and the light dashed contours show the 1, 2 and 3a ranges for the concentration 
distribution. 
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Fig. 2. — The average logarithmic slope rj in the annulus between the images (k oc R 1 ^) for 
adiabatically compressed models of HE2149-2745 as a function of the halo concentration c = r v i T /r s 
and the mass fraction in cold baryons of the visible galaxy /& = Mg /M v i r . The solid lines are 
contours of rj spaced at intervals of At? = 0.2. An isothermal model has rj = 2 (the heavy solid 
contour). The dashed contours show the likelihood of the concentration. The heavy dashed contour 
is the average concentration given the virial mass and the light dashed contours show the 1, 2 and 
3cr ranges for the concentration (see Fig. ||). 
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Fig. 3. — The weak lensing surface density AE at R w = 75h kpc for HE2149-2745 as a 
function of the halo concentration c = r v i r /r s and the mass fraction in cold baryons of the visible 
galaxy fb = Mh /M V i r . The solid lines are contours of AS spaced at intervals of 20hM@/pc 2 
with a heavy contour for the surface density expected from the SDSS weak lensing results, 
AT,sdss — lOO/iM /pc 2 . An isothermal model would have AT,sdss — 7O/iM /pc 2 . Based on 
the virial mass of each model, the dashed contours show the likelihood of the concentration. The 
heavy dashed contour is the average concentration given the virial mass and the light dashed 
contours show the 1, 2 and 3a ranges for the concentration (see Fig. |l]). 
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We consider only the innermost R w = 75h~ 1 kpc bin of the SDSS measurements in order to 
minimize contamination of the signal from sources other than the galaxy halo. In addition to any 
intrinsic uncertainties in the SDSS results, the conversion of the results into a constraint on our 
models is complicated. We assume that we can estimate the central stellar velocity dispersion of 
the lens galaxy using an SIS lens model, which appears to be a fairly reliable assumption (see 
Kochanek et al. |200C| , Treu & Koopmans j2002p . We then use the observed properties of early-type 



galaxies in the SDSS (Bernardi et al. |2002| ) to convert the velocity dispersion into estimates for the 
r' and i' magnitudes the lens galaxy would have locally. This avoids any problems of luminosity 
evolution at the price of introducing the roughly 1 mag of scatter in the Faber-Jackson relation. 



Combining these magnitudes with the r' and i' magnitude zero points from Blanton et al. (2001) 



and a fit to the A£(i?) measurements in McKay et al. ( 2001| ), we estimate that the SDSS weak 



lensing results correspond to a constraint on the aperture mass of 

/ a \3.2 

AZ SDSS (R W = 75/*-^) =i 80 [ 2Q0km/s J hM Q /pc 2 (16) 

with a crudely estimated uncertainty of a factor of two. That the weak lensing signal is dominated 
by the galaxies has been checked by showing that the dynamics of satellites around the target 
galaxies confirm the mass estimates (McKay et al. 2002] ). Future analyses of weak lensing in the 



SDSS, particularly segregating the galaxies by velocity dispersion rather than luminosity, will 
clarify the scalings and the uncertainties. As a reality check of this scaling, a singular isothermal 
sphere of velocity dispersion a predicts AS = 62(<r/2OOkm/s) 2 /iM0/pc 2 , and suggests that the 
scaling of the SDSS results is somewhat high. For now, the weak lensing results provide an 
interesting, but not a definitive, constraint on the halo properties. 



3. Results 



We use the four simple, well-characterized time delay lenses PG1115+080 (Schechter et 
al. [19971 , Barkana gggfl Impey et al. |l998|) , SBS1520+530 (Burud p002b| , Faure et al. pOO^) , 
B1600+434 (Burud et al. p00C| , Koopmans et al. gOOCD and HE2149-2745 (Burud et al. pOOp for 
our analysis, using the Hubble constant scalings derived in Kochanek ( |2002b| ). Of the remaining 
5 time delay lenses, B0218+357 lacks an accurate measurement of the position of the lens galaxy 
(Lehar et al. 2000| ), B1608+656 has two interacting lens galaxies (Koopmans & Fassnacht 1999] ), 
Q0957+561 and RXJ0911+0551 are too heavily perturbed by their parent clusters for our simple 



scaling with the surface density (see Kochanek 2002b ), and the properties of the lens galaxy in 
PKS1830-211 are in dispute (Winn et al. |2002| , Courbin et al. |2002| ). For each lens we model the 
baryonic distribution using the photometric profiles scaled by their intermediate axis scale lengths 
(the geometric mean of the major and minor axes). Since the masses are determined by the lens 
geometry, this is the only observational property of the lens galaxy needed for our calculation. 
The results are insensitive even to large (50%) changes in the scale length. The only important 
uncertainties in the measurements are the time delays and the B1600+434 lens galaxy position. 
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Figures || and |2| show the expected values for (k) and rj in HE2149-2745 as a function of the 
cold baryon fraction fa and the halo concentration c = r v i r /r s . For each value of fb and c, the 
virial mass (and radius) are determined by the observed critical radius of the lens. As the cold 
baryon fraction fa — > 1, the surface density converges to that of a constant M/L model and the 
density exponent is relatively steep (rj ~ 2.5). The images lie somewhat outside the effective radius 
of the lens (R e = (X'48, R c = C/85, Ri = 0'.'34 and R2 = 1'.'35) so the average slope of the surface 
density is closer to the outer density exponent (rj = 4) than the inner density exponent (rj = 1) 
of the Hernquist profile (Eqn. [l]). As we lower the cold baryon fraction fi,, the expected surface 
density (k) rises and the slope rj flattens. They reach the isothermal values, (n) = 1/2 and 77 = 2, 
for baryon fractions fb ~ 0.025. Very extended halos (small c) are too diffuse, and very compact 
halos (large c) are too compact to contribute to (k). 

In simulations, halos with a fixed virial mass are found to have a limited range of 
concentrations (see Bullock et al. 2001 and Eqn. pi]), so only a portion of the fb~c plane can be 



occupied by real halos. We find that the theoretically expected range of concentrations for the 
lenses (c ~ 10) depends little on the cold baryon fraction, as we can see in Figs, [l] and [2] from the 
superposed contours for the expected range of concentrations given by the Bullock et al. ( p001| ) 
results. As a result, the lensing properties depend almost exclusively on the cold baryon fraction. 
If we do not include the adiabatic compression of the halo the results change slightly, but the 
results are similar because the estimates with and without compression must be the same in the 
limits of both large and small baryon fractions. 

Figure || shows the expected strength of the weak lensing signal AS measured R w = 75/i -1 kpc 
from the lens. The aperture mass/surface density which determines the measured weak lensing 
shear steadily declines as we increase the cold baryon fraction and the total halo mass decreases. 
This means that the amplitudes of the shears measured by weak lensing provide a convenient 
observational constraint on the cold baryon fraction. Implementing this in practice is difficult, but 
for the scalings leading to Eqn. (|l|) we predict that AY, S dss - 100/iM Q /pc 2 for HE2149-2745 
given the 214 km/s velocity dispersion estimated from the image separations. With our rough 
estimate of (logarithmic) uncertainties of a factor of two, the discriminatory power of the weak 
lensing estimates is limited but favors models with substantial amounts of dark matter. Like the 
values for (k) and rj, we can largely describe the weak lensing amplitude simply by the cold baryon 
fraction. 

The restricted range of permitted concentrations and the weak dependence of the lens 
properties on the concentration in that range allows us to simplify the results by considering only 
the concentration probability averaged values for (k) and rj as a function of Figs. ||] and || 
show the expected values of (k) and rj as a function of the cold baryon fraction for PG1115+080, 
SBS1520+530 B1600+434 and HE2149-2745. For a fixed cold baryon fraction f b , the results for 
the four lenses are similar, the differences between adiabatically compressed and uncompressed 
models are small and the scatter due to the permitted range of concentrations is still smaller. Over 
the range of baryon fractions from fb = 10~ 3 to = 1, the expected surface density decreases 
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Fig. 4. — The concentration-averaged estimates for 1 — (k) as a function of for PG1115+080, 
SBS1520+530, B1600+434 and HE2149-27145. The solid (dashed) curves are for models with 
(without) adiabatic compression. A horizontal line shows the value expected for isothermal models. 
The diagonally cross-hatched region shows the lower bound on fb based on the local inventory of 
cold baryons by Fukugita et al. ( 1998 ) for Hq = 60 km/s Mpc and S7o = 0.3. The horizontally 
cross-hatched region shows the upper bound on fb set by the global ratio of Qb/fio estimated from 
the CMB. 
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Fig. 5. — The concentration-averaged estimates for the annular density exponent rj (k oc R 1 7? ) 
as a function of f b for PG1115+080, SBS1520+530, B1600+434 and HE2149-27145. The solid 
(dashed) curves are for models with (without) adiabatic compression. A horizontal line shows the 
value expected for isothermal models. The limits on /& are the same as in Fig. 
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Fig. 6. — The concentration-averaged estimates for the ratio AY,/AY<sdss between the expected 
weak lensing signal AS and our best estimate of the SDSS measurements AT,sdss for galaxies 
with the velocity dispersion of the lens. The normalization of AT,sdss is uncertain to a factor of 
2 (0.3 dex). The solid (dashed) curves are for models with (without) adiabatic compression. The 
limits on /& are the same as in Fig. ^j. 
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Fig. 7.— The Hubble constant as a function of f b for the lenses PG1115+080, SBS1520+530, 
B1600+434 and HE2149-27145. The heavy curves show the results for the individual lenses 
including adiabatic compression. The shaded envelope bracketing the curves is the 95% confidence 
region for the combined lens sample. The vertical bands show the lower bound on ff, from 
local inventories, including its H^ 1 scaling, and the upper bound from the CMB. The horizontal 
band shows the estimate of Hq = 72 ± 8 km/s Mpc by the HST Key Project (Freedman et 
al. [2001 ). Without adiabatic compression, the estimates from the time delays predict values for Hq 
approximately 5 km/s Mpc lower. 
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from (k) ~ 0.7 to (k) ~ 0.2 and the exponent of the density profile steepens from r\ ~ 1.6 to 
i] ~ 2.6. B1600+434 behaves differently because it has a smaller ratio between the Einstein radius 
and the half-light radius than the other three systems. 

The cold baryon fraction, /& = M#/M„j r = ^b jC oid/^M, is set by the cosmological density 
^b,coid of baryons that have cooled enough to be modeled by the visible galaxy compared 
to the total density of baryons and dark matter Qm- Local accountings for cold baryons 
(stars, remnants, cold gas components) by Fukugita, Hogan & Peebles (|1998 ) estimated that 



0.0024//i7o ^b,coid ^ O.OO64//170, which sets a Hubble constant-dependent lower bound that 
fb ^ 0.01 for Slo = 0.3. The global ratio of rj^/Oo sets an upper bound on the cold baryon fraction. 
If clusters of galaxies contain a fair sample of material (White et al. 1993j ), their baryonic mass 



fraction determines ^/Qq ~ (0.113 ± 0.005)(1 + 0.22/i 7 q )h 70 ' where the numerical estimate 
is from the recent study of Allen, Schmidt & Fabian ( |2002| ), Current analyses of the cosmic 



microwave background (CMB) anisotropics find Ofe/Oo — 0.15 ± 0.05 (e.g. Netterfield et al. 2002 
Wang, Tegmark & Zaldarriaga |2002| ), consistent with the cluster inventory results. In this context 
we should also note that the CMB analyses generally favor lower Hubble constants than the Key 
Project estimates. The local baryon inventory and the global ratio f^/r^o set upper and lower 
bounds on the cold baryon fraction appropriate for our models. The cold baryon fraction 
could be substantially higher than observed locally at the price of introducing a cold but locally 
unobserved baryonic population. We have superposed these local and global limits on on Figs. |] 
and H 

For baryon fractions similar to local inventories (/;, ~ 0.02), the surface density (k) and 
logarithmic slope r\ are remarkably close to the isothermal values of (n) = 1/2 and 77 = 2. For 
PG1115+080, SBS1520+530, B1600+434 and HE2149-2745 we find («) = 0.50 ±0.05, 0.52 ±0.04, 
0.50±0.03, and 0.52±0.04 respectively, and 77 = 1.98±0.09, 1.95±0.09, 1.91 ±0.05, and 1.95±0.09 
respectively if we constrain ft, by the local values scaled to Hq = 60 km/s Mpc. This means that 
our standard isothermal lens model, which represent our best independent estimate of the mass 
distribution on these scales, is also the expectation from standard CDM halo models. If we allow 
larger cold baryonic mass fractions, then the surface densities decrease and the slopes steepen. We 
find (as) = 0.30 ± 0.04, 0.34 ± 0.03, 0.45 ± 0.11 and 0.35 ± 0.03, and 77 = 2.40 ± 0.06, 2.32 ± 0.07, 
2.02 ± 0.27 and 2.30 ± 0.05 for PG1115±080, SBS1520±530, B1600±434 and HE2149-2745 
respectively when we force fb = f^/f^o = 0.15 ± 0.05. These parameter values lie roughly midway 
between isothermal and constant M/L mass models. 

The expected weak lensing signals, shown in Fig. |6| as the ratio AY,/AT>sdss between 
our estimate from the models and our estimate for the SDSS measurements on a scale of 
R w = 75/i -1 kpc, provides one independent test of the models. The weak lensing signal diminishes 
as fb increases and the total mass of the halo diminishes. Models with f^ similar to local inventories 
predict weak lensing signals roughly consistent with the SDSS measurements (about 70% of the 
expected amplitude, but we only believe the scaling to a factor of two at present). If all available 
baryons cool, the expected weak lensing signal is too weak to agree with the SDSS measurements 
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(only 30% of the expected amplitude). Models with constant M/L predict weak lensing signals 
only 10% of the expected amplitude and are strongly ruled out. With a factor of two uncertainty, 
the weak lensing constraint formally restricts the models to the range 0.002 <S ft, ^ 0.015 (90% 
confidence). Note that when Keeton (2001) used these lens models to fit the distribution of image 
separations in lens surveys, he found it was difficult to match the observed separation distributions 
for f b £ 0.1. 

Finally, in Fig. we use the estimates of (k) and rj for the CDM halo models to estimate 
the Hubble constant as a function of f b for PG1115+080, SBS1520+530, B1600+434 and 
HE1149-2745.Q As we would expect from the homogeneity of the lenses (see Kochanek 2002b| ), 
these four lenses predict similar values of Hq as a function of /;,. For low /&, B1600+434 gives 
somewhat higher values for Hq than the other three lenses. The joint estimate has small statistical 
uncertainties (5%) compared to the overall variable range. If the baryon density is restricted to 
agree with the local baryon inventory, then Hq = 52 ± 6 km/s Mpc (formally at 95% confidence), 
and if it is restricted to agree with the global baryon inventory, then Hq = 65 ± 6 km/s Mpc. 
Restricted by the SDSS weak lensing results rather than f b , we find that Hq = 43 ± 7 km/s Mpc, 
but much of the lower range has baryon fractions below that found in local inventories. If we 
combine the weak lensing constraints with the lower bounds on f b from the local inventories, 
we find that Hq = 48 ± 5 km/s Mpc. The limits for models without adiabatic compression are 
approximately 5 km/s Mpc lower. 

Superposed on Fig. |7| is the local estimate of Hq = 72 ± 8 km/s Mpc by the HST Key Project 
(Freedman et al |2001| ) . The local estimate is consistent with these four lenses only in the limit of a 
constant M/L, as we had found previously (Kochanek [2002ap . Formally, agreement with the Key 
Project requires > 0.06 (one-sided, 95% confidence). It is slightly inconsistent with models in 
which all baryons cool, and grossly inconsistent with models having realistic baryon fractions or 
weak lensing signals consistent with the SDSS measurements. 



4. Discussion 

Because gravitational lens time delays are determined by the Hubble constant and the average 
surface density (k) of the lens galaxy in the annulus between the images (Kochanek p002b| ), we can 
make unambiguous estimates for the behavior of time delays in standard CDM halo models. In 



The scaling solutions from Kochanek (2002b) for these lenses are Ho = A(l — (k)) + B{K,){rj 



1) where the 

coefficients (A, B) in units of km s/Mpc for the lenses are (92.3, 4.6), (93.2, 10.5), (103.9, 20.6) and (84.4, 13.6) for 
PG1115+080, SBS1520+530, B1600+434 and HE1149-2745 respectively, and they produce fractional errors in H 
due to uncertainties in the astrometry and the time delays of 9%, 6%, 9% and 13%, respectively, after broadening 
the errors in the time delays to a minimum of 5% to encompass systematic effects such as convergence fluctuations 



from large scale structure (e.g. Seljak 1994, Barkana 1996). PG1115+080 is in a group which supplies an additional 
external convergence of K ext = 0.2(1 — (k)). 



-17- 



these models, the expected delay is controlled by the mass fraction, fb = Mn/M v i r = Qb,cold/&o, 
in cold baryons making up the observed lens galaxy relative to the overall halo. As the cold baryon 
fraction rises, so does the Hubble constant. 

When the cold baryon fraction is comparable to the local baryonic content of galaxies 
(fb — 0.02, Fukugita et al. 1998| ), the model parameters closely match those for isothermal (flat 



rotation curve) dark matter dominated lens models and the halos produce weak lensing signals 
compatible with weak lensing measurements in the SDSS (McKay et al. [200 1| ). The mean surface 
density in the annulus is almost exactly («) = 1/2 and the local slope of the surface density is 
almost exactly rj = 2 (k oc R 1 ^ 11 ). Isothermal models are not only the best observational estimate 
for the lensing potential on these scales, they are also the model predicted by CDM assuming 
standard parameters and baryonic populations. For baryon fractions with a lower limit set by 
the local inventory and the upper limit set by the weak lensing measurements, we find that 
H = 48 ± 5 km/s Mpc based on the time delays measured for PG1115+080, SBS1520+530, 
B1600+434 and HE2149-2745. If all baryons were to cool and f b ~ 0.15 ± 0.05, based on 
constraints from either the CMB (e.g. Netterfield et al. 2002, Wang et al. 2002j ) or cluster baryon 



fractions (e.g. White et al. 1993| , Allen et al. 2002j ), then the Hubble constant could be as high as 



Hq = 65 ± 6 km/s Mpc. Such models require most of the cold baryons in the lens galaxies to be 
in a locally invisible population and correspond to mass distributions less consistent with direct 
estimates. Both of these possibilities are lower than the local estimates of Hq = 72 ± 8 km/s Mpc 



from the HST Key Project (Freedman et al. 2001 ), which agrees with the time delays of these 



four lenses only for mass distributions with constant M/L ratios. Thus, our detailed models for 
the expected properties of time delays in standard CDM halos agree with our simple models in 



Kochanek (2002a), and we are faced with a conflict between CDM halo models, gravitational lens 



time delays and the local distance scale. 

While there is some room for error in the lens results, the mutual agreement of the four simple, 
well-characterized time delay lenses and the simple relation between time delays, surface densities 
and the Hubble constant makes it difficult to point to a weakness (Kochanek |2002b ). The most 



important observational steps are to improve the accuracies of the existing delay measurements 
and to expand the number of systems with delay measurements. If the homogeneity of the results 
for simple lens systems, as compared to more complicated systems in clusters or with interacting 
galaxies, continues, the case for the existence of a conflict will become overwhelming. Improved 
characterizations of the lenses, either to constrain the mass distribution in the time delay lenses 
directly or to allow us to include the five other time delay lenses, are also important, but depend 
on obtaining deeper HST imaging of the systems. Other constraints on the mass distributions 
such as weak lensing or the stellar dynamical measurements of the lens galaxies can also help to 
break any degeneracies. In particular, estimates of the weak lensing signal as a function of the 
stellar velocity dispersion rather than luminosity would be excellent constraints on the halo extent 
in time delay lenses. 



The systematic uncertainties in the mass distribution can be minimized by measuring time 
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delays in lenses where the baryons dominate the mass and there there is little difference between 
a constant M/L model and a model with dark matter. This means measuring the time delay 
in very low redshift lens galaxies where the ratio of the critical radius to the effective radius, 
R c /R e oc Dol is small and the mass near the Einstein ring is increasingly dominated by the 
baryons. For example, models of Q2237+0305 at z\ = 0.04 suggest that less than 10% of the 
mass inside the Einstein ring of the lens can be dark, instead of the roughly 50% for typical 
models of higher redshift lenses (Trott & Webster |2002| ). Unfortunately, Q2337+0305 has shown 
no variability on the very short time scale of its expected delay, making it a poor candidate for 
measuring time delays. There is a certain irony to proposing that local galaxies, which might be 
incorporated in local distance scale studies, are the ideal time delay lenses, but it may also lead to 
a system where the local and the "cosmological" distance scales can be compared directly. 

Acknowledgments. CSK thanks D. Rusin, P. Schechter, U. Seljak, J. Winn and S. Wyithe for 
discussions and comments. CSK is supported by the Smithsonian Institution and NASA ATP 
grant NAG5-9265. 



A. Tidal Truncation 

Many lens galaxies, including many of the time delay lenses, are group and cluster members 
(e.g. Keeton, Christlein & Zabludoff 200C| ). Particularly when PG1115+080 appeared to be 



the exception among time delay lenses rather than the rule, tidal truncation of its halo by the 
surrounding group was frequently invoked as a possible explanation for what appeared to be an 
anomalously low Hq estimate for isothermal lens models (e.g. Impey et al. 1998| , Koopmans &; 



Fassnacht |1999| ). The idea is that a tidally truncated halo is more centrally concentrated and 
has a lower (re), so it produces a higher estimate of Hq for a fixed time delay. However, in order 
to significantly raise Hq, the halo must have a tidal truncation radius comparable to the critical 
radius of the lens and the surface density of the truncating halo must be kept low. While orbital 
history, phase, projection effects and differences in the density distribution make the physics of 
tidal truncation considerably more complicated then the simple model we discuss, simple scaling 
laws for tidal truncation suggest that this standard picture is physically implausible. 

The first problem with tidal truncation is that it is extraordinarily difficult to truncate the 
mass distribution of a lens galaxy on scales comparable to the critical radius. For two SIS lenses, 
the tidal radius of the less massive lens is Rt = Rp^o/bp) 1 / 2 when the projected separation is 
equal to the true separation in three dimensions. The ratio roughly equals the velocity dispersion 
ratio of the two halos, (bo/bp) 1 / 2 ~ 0o/<7 p , so for the typical massive lens galaxy in a low mass 
group or cluster the tidal radius Rt <^ R p /2. Truncating the lens galaxy near the Einstein ring 
is almost impossible because the inequality Rt < &o implies a galaxy inside the critical radius of 
the cluster, R p < {b p bo) l l 2 < b p , where it cannot produce the image geometries seen for isolated 
lenses. In fact, weak lensing studies of galaxies in the cores of rich clusters have shown that even 
in these high density environments the halos of the galaxies seem to be truncated only on scales 
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J> 20/i 1 kpc far larger than the ^ 5h 1 kpc scales on which tidal begins to alter the time delays 



significantly (Natarajan et al. 2001). 



The second problem is that the reduction in (k) from tidally truncating the primary lens 
must be accompanied by an increase in (k) from the convergence K ex t of the halo responsible 
for the truncation. Existing lens models partly incorporate this effect through the scaling of 
external shears with surface density. As a concrete example, suppose we have a cruciform lens 
with images on the symmetry axes. We consider two models: an SIS primary lens perturbed by 
an SIS cluster, and a point mass primary lens perturbed by an SIS cluster. We model the cluster 
as an equal external shear and convergence, 7 = n = b p /2R p , for a cluster with critical radius 
bp at impact parameter R p . For a source directly behind the lens, the two models can produce 
identical image positions. The critical radius of the primary lens, bo, is the same in both models, 
but the cluster shear for a isothermal primary is less than that for a point mass primary, with 
Ipnt = i^pnt = 2ksjs(1 — Ksis) = ^Isisi^ — Jsis)- As the lens galaxy becomes more compact, the 
external shear and convergence from the cluster increase. The time delay ratio depends only on 
the strength of the cluster perturbation, 



At pn t _ 1 — KSIS j n 
AtsiS Ksis 



1 



1 - 2ks/5 



2(1 - ksis) + ■■■ (Al) 



for small perturbations. The factor of 2 is the ratio of the (k) values for a point mass and an 
isothermal primary lens, and the correction is due to the higher convergence of the cluster when 
the primary lens is a point mass. The delay ratio is monotonically declining, reaching unity for 
K sis = ISIS = 0.36 ( r y pn t = 0.45). Thus, as we raise the cluster convergence, so as to make 
tidal truncation physically plausible, the reduction in Hq by the cluster convergence eventually 
becomes more important than the increases in Hq from making the primary lens more centrally 
concentrated. 

While the rise in K ex t produced by the natural scaling of lens parameters with the increasing 
concentration of the lens galaxy mimics the effects of tidal truncation, it does not necessarily 
lead to self-consistent results because the amplitudes of the shear and convergence, j ex t and K ex t, 
from the perturbing halo are closely related to the tidal truncation radius. For an isothermal 
galaxy truncated by an isothermal halo, the surface density of the perturber at the tidal radius is 
Kext = bo/2R t — > 1/2 as the tidal radius approaches the critical radius, Rt — > 60 • Equivalently, 
there is a close relation between the tidal radius and the external shear from the perturbing halo, 
with Rt = (bobp) 1 / 2 /2j ex t for two isothermal halos, which means that the truncation radius cannot 
be close to the Einstein ring unless the external shear is large. Hence, the standard parameters of 
"tidally truncated" models of PG1115+080, where K ex t = jext — 0.15 are inconsistent with simple 
models of the group parameters required to tidally truncate the lens galaxy. The standard model 
will tend to overestimate Hq for tidally truncated galaxies by underestimating K ex t- 
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